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Abstract 

Explicit expressions for one point moments corresponding to stochas- 
tic Verhulst model driven by Markovian coloured dichotomous noise 
are presented. It is shown that the moments are the given functions 
of a decreasing exponent. The asymptotic behavior (for large time) of 
the moments is described by a single decreasing exponent. 

Keywords: Stochastic Verhulst model, one point moments, explicit 
expressions. 

1 Introduction 



There are a lot of papers devoted to description of the temporary evolu- 
tion of moments with exactly solvable nonlinear stochastic equations. In [1] 
we gave some general procedure to explicitly solve the master equations of 
hyperbolic type corresponding to nonlinear stochastic equations driven by 
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dichotomous noise. The method is based on a generahzation of Laplace fac- 
torization method [21 E]- As an example we have considered complete exact 
nonstationary solution of the master equations for probability distribution 
corresponding to stochastic Verhulst model. 

In this paper we calculate one points moments of arbitrary degree and 
discuss its time evolution. Let us consider nonlinear stochastic dynamical 
system 

X = p{x) + a{t)q{x), (1) 

where x{t) is the dynamical variable, p{x), q{x) are given functions of x, a{t) 
is the random function with known statistical characteristics. The model 
([1]) arises in different applications (see for example [H [5] and bibliography 
therein). An important application of this model consists in study of noise- 
induced transitions in physics, chemistry and biology. The functions p{x), 
q{x) are often taken polynomial. For example, if we set p{x) = pix +p2x'^, 
q{x) = q2X^, pi > 0, p2 < 0, \p2\ > q2 > 0, then the equation ([T]) describes the 
population dynamics when resources (nutrition) fluctuate (Verhulst model). 
In the following we will assume a{t) to be binary (dichotomic) noise a{t) = 
±1 with switching frequency 2i> > 0. As one can show (see [6j), the averages 
W{x,t) = {W{x,t)) and Wi{x,t) = {a(t)W{x,t)) for the probability density 
W{x, t) in the space of possible trajectories x(t) of the dynamical system ([T]) 
satisfy a system (also called "master equations"): 

+ ipix)W)^ + iqix)W,)^ = 0, 
{W^)t + 2uWi + {p{x)W^)^ + {q{x)W)^ = 0. 

We suppose that the initial condition W{x,Q) = Wo{x) for the probability 
distribution is nonrandom. This implies that the initial condition for Wi {x, t) 
at t = is zero: Wi{x,0) = {a{0)W{x,0)) = {a{0))Wo{x) = 0. The prob- 
ability distribution W{x,t) should be nonnegative and normalized for all t: 
W{x, t) > 0, JZo W{x, t) dx = 1. 

In |lj we have obtained the following explicit form of the complete solution 
of the system ([2]) for probability distribution W{x,t): 

W{x, r) = ie-^ {S {x - j^^^f^^^^) + S {x - } + 

2^ {e^il+{p2-q2)-x{p2-q2)) - ^ H ( e-(l+(p2+g2)-2:(P2+92)) ^ ^*) } ' 

(3) 

where t = vt is the dimensionless time and x^, is an initial value for ([T]), 
H{z) = /f 5{6)d0 is the Heaviside function. 

The solution ([3]) corresponds to Cauchy problem Wo{x) = 5(a; — x*). Here 
we set that u = 1 . The function W{x,t) is in fact a conditional probability 
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distribution, that is W{x,t)Ax = W{x,t \ x^:,t = 0)Ax is the probabihty 
that at the time r the dynamical variable x belongs to interval {x, x + Ax) 
under condition that at some previous initial time r = the variable x is 
equal to x*. 

From the equation ([T]) follows that the dynamical variable has three sta- 
tionary points: 

1 1 

Xi = — , X2 = . , X3 = 0. 

\P2\ + q2 \p2\-q2 

It is convenient to use the definition xi and X2 for transformation of the 
expression ([3]) to the form: 

Wix, t) = -e <6 \x z ^ — ] + 6 \ x z ^ — } + 



X1X2 / xxie \ TT ^^2e 

H : TT - x^ \ - H \ — - X* 



{x2 — xi)x^ [ \xi + x{e "^ — 1) / \X2 + x{e ^ — 1^ 

2 Calculation of one point moments 

The one point conditional moments of n-th order one defines as 



(4) 



kJt) = (x"(r)|x(0) = x*,r = 0) = / x"l^(x,r)rfx, (5) 

J{D) 

where (D) is the support of the probability distribution. Further we consider 
the case (D) = {xi,X2). After simple calculations one obtains 

-L J- I / 6 X^X2 \ I 6 X^X^ 



""-^^^ - 2''~\[x2 + (er - l)xj +l^xi + (e--l)xj / + 



(6) 



- — — - {{X2mr^' - {xMr)r-'} 

[X2 - Xij{n - 1) 



where 

Pir) 



7(r) 



X* + {x2 - x*)e" 



3 



Let r ^ 0, then /3(t) — > ^ and 7(r) — > In this hmit from one has 
i^n{T) ^l- Let us consider another asymptotic r ^ oo. From ^ one 
obtains for n = 1 the following stationary value of the moment 

Ki(r) = ^ (lnx2 — Inxi), 
and for n ^ 1 stationary values of moments are equal to 



n 



Generally the moments Kn(r) are given functions depending on the de- 
creasing exponent e~'^ and can be represented by a series over the powers 
of e""^. In [H] a similar representation was found for the stochastic Verhulst 
model with fluctuating coefficient at the first degree of the dynamical variable 
X. In the limit for large r ^ 1 the time behavior of K„(r) can be described 
by a single exponent. 

Important role is played by the first two initial moments (n = 1,2). Let 
us consider the moment of first order. In this case 

X1X2 X2I3{t) 




X2-X1 Xi7(r)' 

In the asymptotics r ^ 00 from ([7]) in first order over the infinitesimal 
exp(— r), one has 

, , X1X2 , X2 / Xi -\- X2 X1X2' 

Ki{t) ^ In h 



X2 — X\ X\ \ 2 x^ / 

It is interesting that there exists the initial value = = In 

this case the coefficient at exp(— r) is equal to zero. Therefore we should 
take into account the next order, i.e. exp(— 2r). Physically it means that in 
point x^, = the correlation of variable x(r) with the given initial value of 

variable x (x(0) = x*) decreases more rapidly at r ^ 00. When x* the 
correlations tends to stationary level as exp(— r). 

Here we give also an explicit expression for the case n = 2: 



«2(r) = ^e" 



^X2 + (e^ - l)x*/ \xi + (e^ - l)x^ 

x^XiX2(l - e~^) 

(x* + (x2 - x*)e-^)(x* - (x* - xi)e"^) ' 



(9) 



4 



3 Concluding remarks 



We have considered the time evolution of one point moments of dynamical 
variable corresponding to the stochastic Verhulst model. The explicit form of 
the moments shows that the moments are the functions of exp(— r). In [8j it 
was shown for Verhulst model when parameter fluctuates at dynamical vari- 
able X (not x^), that the exact solution for one point moments is presented 
by a series over powers of exp(— r). From formulae obtained here one can 
write the moments Knir) in the same form. It should be noted that in this 
communication we have obtained an explicit form of the solution. Under 
the condition r — oo the moments decrease in time as single exponential 
function. It is shown that the time dependence of the moment k,i{t) which 
physically describes the correlation between the value of the dynamical vari- 
able X at the time r with its given (nonrandom value x^) at the initial time 
r = changes. This time behavior depends on the choice of the initial value 
x^,. The critical value is = l/|p2|- For this value of x^, in the limit t ^ oo 
the correlations decrease as exp(— 2r), not as exp(— r). 

It should be noticed that the one-point moments for some special type of 
the dynamical system ([T]) with polynomial functions p{x) and q{x) Gaussian 
white noise fluctuation coefficient at the first power of x were considered in 
[71 El HBJ , where it was shown that the asymptotic behavior of the moments 
is described by a power function. 
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